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THE ORIGIN OF BLUE CORES IN HUBBLE DEEP FIELDS E/SO GALAXIES. 
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ABSTRACT 

In this letter we address the problem of the origin of blue cores and inverse color gradients in early- 
type galaxies reported in the Hubble Deep Field North and South (HDFs) by Menanteau, Abraham, 
& Ellis (2001). We use a multi-zone single collapse model. This model accounts for the observed blue 
cores by adopting a broad spread in formation redshifts for ellipticals, allowing some of these galaxies 
begin forming no more than ~ 1 Gyr before the redshift of observation. The single-zone collapse model 
then produces cores that are bluer than the outer regions because of the increase of the local potential 
well toward the center which makes star-formation more extended in the central region of the galaxy. 
We compare the predicted Vqqq — isuM color gradients with the observed ones using the redshift of 
formation (zp) of the elliptical as the only free parameter. We find that the model can account with 
relatively good agreement for the blue cores and inverse color gradients found in many spheroidals and at 
the same time for the red and smooth colors profiles reported. Based on the model our analysis suggests 
two populations of field ellipticals, one formed recently, within < lGyr and another much older formed 
> 4Gyr since the redshift of observation. 

Subject headings: galaxies: elliptical, lenticulars — galaxies: formation — galaxies: evolution 



1. INTRODUCTION 

The study of field elliptical galaxies has seen a dramatic 
boost in the past couple of years, mainly from the real- 
ization that these systems can be objectively used to test 
opposed models of galaxy formation. 

For a number of years it was accepted the hypothesis 
that field elliptical galaxies formed in isolation during a 
single burst of formation at high redshift. This had its ori- 
gin in the early monolithic collapse model postulated by 
Eggen et al. (1962) and later sustained by observations of 
ellipticals in rich clusters over a wide redshift range which 
closely follow a number of "fundamental relations". The 
low scatter in the fundamental plane (Djorgovski & Davis 
1987; Dressier et al. 1987) and in the color-magnitude re- 
lation (CMR) (Sandage & Visvanathan 1978; Bower et al. 
1992; Ellis et al. 1997; van Dokkum et al. 1998) implied a 
high degree of homogeneity in the stellar population, both 
consistent with a model in which ellipticals formed at high 
redshift and since then passively evolved. 

On the other hand, in hierarchical CDM models (White 
& Rees 1978; White & Frenk 1991) galaxies are formed 
only as the result of merging of smaller sub-units, suggest- 
ing that massive ellipticals could only have been assembled 
recently (since z ~ 1) and in consequence there would be a 
paucity of early-types at higher redshifts (Kauffmann et al. 
1996). Unfortunately, much of the observations were early 
focused on elliptical galaxies in rich clusters, where the 
existence of old ellipticals at high redshift is expected in 
hierarchical CDM models because evolution is accelerated 
in these dense regions (Governato et al. 1998; Kauffmann 
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et al. 1999). 

Several studies have attempted to test whether field el- 
lipticals were formed at high redshift and evolved in isola- 
tion since, following to some extent the methodology from 
Moustakas et al. (1997) and Zepf (1997), consisting of us- 
ing a combination of optical and infrared colors for select- 
ing spheroidals in search for a decline in the co-moving 
density of spheroidals; all of them with consistent results, 
a significant lack of red objects compared with the predic- 
tions of pure luminosity evolution (PLE) models (Barger 
et al. 1999; Daddi et al. 2000). The addition of HST imag- 
ing lead to the conclusion that a model with constant co- 
moving density and a single epoch of formation for E/SO 
at all epochs is ruled out. A more conclusive test in this 
line will need a complete census of ellipticals up to z = 1, 
recent results from Im et al. (2000) are inconclusive but 
compatible with very mild density evolution. Similarly, 
using dataset multicolor information it has been reported 
a population of E/SO with colors significantly bluer than 
passive evolution (Kodama et al. 1999; Abraham et al. 
1999). In a recent detailed analysis Menanteau, Abraham, 
& Ellis (2001) reported, using resolved multicolor photom- 
etry of E/SO in both HDF- North and South, solid evidence 
for the continuous star- formation in E/SO galaxies. The 
resolved ( V — I) data revealed a population of spheroidals 
with central blue cores, color inhomogeneities and inverse 
color gradients, which were indicative of formation activ- 
ity at the galaxies redshift of observation. Interestingly, a 
subset of spheroidals had the color properties of old stel- 
lar population in good agreement with the ones prescribed 
with a stellar population formed in a single burst at z > 2. 

In recent years, the hierarchical CDM semi-analytical 
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Fig. 1. — An example of the ability of the model to repro- 
duce a range of different color gradients. Figure shows from left 
to right: /8i4~band surface brightness map, Vqoq — Isi4 color pixel 
map and Vqoq — 1814(f) color gradient for several E/SO in the sam- 
ple. Open circles represent observed gradients while solid lines 
the model predictions. The complete set of fits can be found at 
http : //www . ociw . edu/~f elipe/blue-cores 



models have dominated the predictions in the field leaving 
alternative frameworks relatively unexplored in which to 
interpret data. The realization that spheroidals, even at 
intermediate redshifts, are a mixture of both a mild star- 
forming and old stellar population, makes imperative to 
search for alternative formation models which can success- 
fully reproduce both properties, the blue cores and smooth 
color profiles observed. In the past, some authors have in- 
vestigated the effect of adding small amounts of SF activ- 
ity to account for some of the blue colors observed in E/SO 
(e.g. Menanteau et al. 1999, Trager et al. 2000). Although 
they have been successful to some extent in reproducing 
the observed colors, they rely on somewhat ad-hoc recipes. 

Recently, Friaga & Terlevich (2001) have used a multi- 
zone chemo-dynamical model to predict the color gradi- 
ents of spheroids at different redshifts. In their model an 
important component to reproduce the blue cores is a late 
infall of gas that extents the duration of the star formation. 
In this letter we investigate a similar model but without 
the late infall feature and take it a step further by produc- 
ing detailed fits to the observed color gradients of all 77 
E/SO studied by Menanteau, Abraham, & Ellis (2001) in 
order to asses whether such single-collapse models really 
fit the observations at all radii and not only the qualitative 



trend of a blue core. We show that the color gradients are 
determined by the shape of the gravitational potential and 
therefore the different escape velocities of the gas rather 
than the need of late infall. In our model the majority of 
the gas (> 80%) is assembled at a very early epoch in the 
formation of the elliptical. This model could be represen- 
tative of how ellipticals assemble in a warm dark matter 
scenario (where small scale structure is suppressed) or in 
a CDM scenario in which the subsequent merging of small 
objects on a big galaxy has been suppressed -maybe due 
to the lack of small scale power (e.g. self-interacting dark 
matter) or strong feedback from the big central galaxy. 

2. THE HDF NORTH AND SOUTH E/s0 SAMPLE 

The E/SO used in our analysis are taken from the field 
sample of galaxies selected by Menanteau, Abraham, & 
Ellis (2001) from the HDFs. Here we include a brief de- 
scription of it, but we point the reader to Menanteau, 
Abraham, & Ellis (2001) for a complete description of it. 
The sample consists of 77 spheroidal galaxies, classified 
as E, E/SO and SO under the Hubble scheme with a limit 
near-total magnitude of /gi4 = 24 mag. Throughout this 
letter magnitudes are in the Vega system. Source detec- 
tion and integrated photometry were taken from Sextrac- 
tor catalogs constructed by the Space Telescope Science 
Institute following the IAU IDs given there. The morpho- 
logical selection of E/SO was performed using both visual 
and machine based classification in order to ensure a ro- 
bust morphology selection (see Menanteau, Abraham, & 
Ellis 2001). Redshift information was taken from publicly 
available compilations, and augmented with photometric 
estimates provided by S. Gwyn (private communication). 

2.1. Observed color gradients 

The quantity that we use to compare the properties of 
spheroidals with the model predictions (see next section) is 
the color gradient of galaxies. We focus only on the Vqoq — 
2"8i4 colors to estimate the color gradient Vqoq — 2su(r) 
as a function of the radius r, given that these represent 
the bands with higher signal-to-noise in both HDFs. In 
order to obtain the color gradients, we first compute the 
centroid and second order moments (x 2 , y 2 , xy) of the 
galaxy utilizing the 2si4-band, from which we estimate the 
ellipticity parameters of the objects (A,B,0). It is worth 
noting that we keep the same elongation ratio e = A/B 
up to the isophotal limit to which the color gradient is 
computed. Finally using concentric ellipses, we compute 
the Vqoq — 28i4 (r^) gradients as the median color in the 
shell between n-\ and r^+i. 

3. MAKING AN ELLIPTICAL 

We model elliptical galaxies as a system with spherical 
symmetry and multiple zones. In particular it is assumed 
that the bulk (> 80%) of the gas in this model was in place 
at the time of formation and was able to form stars, i.e. 
was cool enough. The galaxy is then divided in spherical 
shells -100 for the present case- each of them independent, 
i.e. no transfer of gas is allowed among shells. In each of 
these shells star formation proceeds according to a Schmidt 
law: SFR= vp z ^(t\ where p gas (£) is the volume gas den- 
sity in the shell and v = 8.6(M gas /lO 12 M )-° 115 Gyr- 1 . 
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Fig. 2. — Time since the most likely redshift of formation (zp) 
for all E/S0 in the HDFs sample as a function of its redshift of 
observation, zp is calculated as described in the text. Error bars 
represent la values on the x 2 fitting. 



The initial mass function is assumed to be a power law 
{4>{m) oc m -0,95 ). Star formation proceeds in each shell 
until the gas is heated up by SN to a temperature T that 
corresponds to the escape velocity of each shell (see Mar- 
tinelli, Matteucci, & Colafrancesco 2000 for a detailed de- 
scription of the model). The gas in the elliptical galaxy is 
assumed to be within a dark matter halo of mass 7 times 
larger than the gas mass (Q m = 0.35 and ^5 = 0.05). The 
dark matter follows the density profile described in Mar- 
tinelli, Matteucci, & Colafrancesco (2000). The chemical 
enrichment of the gas and stars is followed in detailed us- 
ing the stellar yields of Woosley (1987) for massive stars 
and supernovae II (M > 8M ), Renzini & Voli (1981) for 
low and intermediate mass stars (0.8 < M/Mq < 8) and 
Thielemann, Nomoto, & Hashimoto (1993) for SN type la; 
we also took into account the lifetimes of stars. For each 
shell we assume that mixing of the gas is very efficient in 
the whole shell and shorter than the lifetime of the most 
massive stars. 

For different masses the model predicts a different time 
dependence of the SFR. In more massive systems the po- 
tential well will be deeper and therefore it will take longer 
for the gas to reach temperatures larger than the escape ve- 
locity in the potential, thus star formation will last longer 
than in less massive systems. Also, for a fixed mass, since 
the potential is deeper in the core of the galaxy, the model 
predicts that star formation will last longer in the core 
than in the outer regions (see Martinelli, Matteucci, & Co- 
lafrancesco 2000). The model provides us with the mass 
fraction that is turned into stars at a given time and their 
chemical composition. This output is then fed into a syn- 
thetic stellar population code (Jimenez et al. 1998) to com- 
pute the spectral properties of the galaxy at different radii 
and times. 

4. RESULTS 

The aim of our analysis is to directly compare the ob- 
served 1^06 — ^814 (r) color gradients with the predictions 
of the above model. In the first place we are interested to 
see whether the color gradients can be reproduced at all. 



We note that the only free parameter now in our model 
is the redshift of formation given that mass is determined 
by the brightness of the galaxy at the redshift of obser- 
vation. To accomplish this we use a maximum likelihood 
method (x 2 ) to compute the most-likely redshift of forma- 
tion (z h p st ) for the best-fitting model. In order to avoid 
spurious results from small fluctuations in the gradients 
due to noise, the model and data gradients were re-binned 
to a common grid of 5 — 6 shells up to a maximum phys- 
ical radius of 10 kpc. In order to transform the observed 
color gradients to physical (kpc) length, we assumed a flat 
cosmology with Q m = 0.35 and Ho = 65 km s _1 Mpc -1 . 
Finally, as the HST point spread function (PSF) is wave- 
length dependent, in order to properly compare observed 
and modeled gradients, we convolved each modeled surface 
brightness profile with a ID radial Gaussian PSF of dif- 
ferent FWHM in the ^06 and i8i4~band (Casertano et al. 
2000). 

We applied this methodology to the whole sample of 77 
E/S0 galaxies. To illustrate the method, in Figure 1 we 
show a selection of four representative galaxies in the sam- 
ple. We can see how the model can successfully account 
for the observed range of Vqo6— ^814(0 color gradients. The 
complete set of fits can be found at http : //www . ociw . edu/~f elipe 
It is worth noting the ability of the model to reproduce the 
blue cores (inverse steep gradients) present in the sample 
(Figure 1 upper two panels), both the color difference and 
the physical scale at which this occurs. In addition very 
flat and smooth color profiles can be successfully repro- 
duced (Figure 1 lower two panels) as expected from a stel- 
lar population formed in single burst at high redshift. To 
emphasize the power of our model, we focus on the bottom 
panel elliptical at z = 0.089. This exhibits a slightly red 
color profile induced by a metallicity gradient, which the 
model naturally accounts for due to the different rate at 
which star formation proceeds in each shell yielding to a 
metallicity gradient. 

The redshift of formation is the relevant measurement 
in our analysis. This is of importance to understand when 
E/S0 assembled to their final morphological shape in which 
we see them today. Figure 2 shows the time since z b p st 
(computed galaxy redshift of formation) as a function of 
the observed redshift for the whole sample. The most 
striking feature is the bimodal distribution of times: a 
population seems to be formed within 1 Gyr, while a sec- 
ond population seems to have formed at least 4 Gyr since 
they were observed. We note that for the later population 
the time since formation could in principle be bigger than 
4 Gyr, since the colors at old ages do not contain enough 
information as to give a precise measurement of the red- 
shift of formation. This population, which the model pre- 
dicts was assembled at an early epoch, agrees with the 
passive evolving population detected independently under 
the basis of S(V — I) color scatter by Menanteau, Abra- 
ham, & Ellis (2001). It is worth to notice that the model 
yields a relatively low redshift of formation for some galax- 
ies when they have large inverse color gradients or marked 
blue colors in order to account for them. 

Figure 3 shows the distribution of the formation red- 
shift for the whole sample. The top axis shows the look- 
back-time. About 25% of field E/S0 in our sample have 
formed at z > 4, with ~ 30% of the sample having formed 
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Fig. 3. — Histogram for the most likely redshift of formation, zp 
values for the whole sample of E/SO. 



at z < 1. The medium redshift of formation is z ~ 2 
and therefore the medium age of the field ellipticals in 
the HDFs is ~ 11 Gyr. Therefore, as a whole, field el- 
lipticals are only 1-2 Gyr younger than cluster ellipticals, 
in agreement with the findings by Bernardi et al. (1998) 
who compared the Mg2 — cr relation for field and cluster 
ellipticals. The main feature of Fig. 3 is the continuous 
formation of field E/SO with redshift, we do not find evi- 
dence for a single epoch of spheroid formation. Although 
for the reasons discussed above the redshift of formation 
can be quite uncertain at the high-redshift end, it is clear 
that we can distinguish a trend of a bimodal distribution in 
the formation redshift. This may be indicative of two pre- 
ferred epochs of formation, according to the model. How- 
ever, this result should not be over interpreted, as this is 
a model dependent one based on the success of the model 
to reproduce the observations. 

5. CONCLUSIONS 

In this letter we have presented an alternative approach 
for the formation of E/SO galaxies. Recognizing the fea- 
tures that E/SO show at low and intermediate redshift 
our analysis seeks to exploit the resolved information now 
available for E/SO to quite high redshift focusing on their 
color gradients as the basis of our analysis. In order to 
obtain information from the color-gradients about the for- 
mation epoch of E/SO, we have adopted a multizone single 
collapse model for the formation of spheroids. In the con- 
text of this model, our main conclusion is that E/SO can be 
properly described by a model in which 80% of the gas was 
processed at an early stage, but in which spheroidals are 
formed continuously as a function of redshift. Although 
there is no preferred epoch for their formation, our analy- 
sis suggest the presence of two population, one assembled 
recently, within ~ lGyr and another with the properties 
of an old stellar population formed > 4Gyr from the red- 
shift of observation. It remains to be learned what star 
formation histories can be obtained in a model indepen- 
dent fashion from the spectra of these objects and whether 
this yields to similar conclusions. 
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